﻿clc,clear
close all
disp("intrinsic function system");

x = 0:pi/1000:2*pi;

function y1r = y1(x, n)
    y1r = sin(n*pi* x /1);
end

function y2r = y2(x,n)
    y2r = cos(n*pi*x/1);
end

function y3r = y3(x,n)
    y3r = sin((n+1/2)*pi* x /1);
end

function y4r = y4(x,n)
    y4r = cos((n+1/2)*pi* x /1);
end

figure(1)
%set(gcf,'color',[1 1 1])
set(gcf,'position',[150 10 1000 800])

subplot(221);
plot(x, y1(x,1),'r');hold on;
plot(x, y1(x,2),'g');hold on;
plot(x, y1(x,3),'b');hold on;
plot(x, y1(x,4),'k');hold on;
title(['cos(n*pi*x/1)'], 'Fontsize',15);
legend('n=1','n=2','n=3','n=4');

subplot(222);
plot(x, y2(x,1),'r');hold on;
plot(x, y2(x,2),'g');hold on;
plot(x, y2(x,3),'b');hold on;
plot(x, y2(x,4),'k');hold on;
title(['cos(n*pi*x/1)'], 'Fontsize',15);
legend('n=1','n=2','n=3','n=4');

subplot(223);
plot(x, y3(x,1),'r');hold on;
plot(x, y3(x,2),'g');hold on;
plot(x, y3(x,3),'b');hold on;
plot(x, y3(x,4),'k');hold on;
title(['sin((n+1/2)*pi* x /1)'], 'Fontsize',15);
legend('n=1','n=2','n=3','n=4');

subplot(224);
plot(x, y4(x,1),'r');hold on;
plot(x, y4(x,2),'g');hold on;
plot(x, y4(x,3),'b');hold on;
plot(x, y4(x,4),'k');hold on;
title(['cos((n+1/2)*pi* x /1)'], 'Fontsize',15);
legend('n=1','n=2','n=3','n=4');